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ABSTRACT 

We present VIMOS integral field spectroscopy of the brightest radio-quiet QSO on the southern sky HE 1029-1401 at a redshift 
of z = 0.086. Standard decomposition techniques for broad-band imaging are extended to integral field data in order to deblend the 
QSO and host emission. We perform a tentative analysis of the stellar continuum finding a young stellar population (< lOOMyr) or 
a featureless continuum embedded in an old stellar population (lOGyr) typical for a massive elliptical galaxy. The stellar velocity 
dispersion of cr^ - 320 + 90km/s and the estimated black hole mass log(MBH/^o) - 8.7 + 0.3 are consistent with the local Mbh-c^* 
relation within the errors. For the first time we map the two-dimensional ionised gas distribution and the gas velocity field around 
HE 1029-1401. While the stellar host morphology is purely elliptical we find a highly structured distribution of ionised gas out to 
16kpc from the QSO. The gas is highly ionised solely by the QSO radiation and has a significantly lower metallicity than would be 
expected for the stellar mass of the host, indicating an external origin of the gas most likely due to minor mergers. We find a rotating 
gas disc around the QSO and a dispersion-dominated non-rotating gas component within the central 3 kpc. At larger distances the 
velocity field is heavily disturbed, which could be interpreted as another signature of past minor merger events. Alternatively, the 
arc-like structure seen in the ionised gas might also be indicative of a large-scale expanding bubble, centred on and possibly driven 
by the active nucleus. 

Key words. Galaxies: active - Galaxies: ISM - quasars: emission-lines - quasars: individual: HE 1029-1401 



1. Introduction 

Since the underlying host galaxie s of seve r al qua si- stellar ob- 
■ jects (QSOs) were resolved by iKristianI (Il973h . QSO host 
galaxies have been extensively studied to understand their 
properties and their relevance for the evolution of the over- 
all galaxy population over cosmic time. Most of our knowl- 
\ edge of QSO host galaxies comes from numerous ground and 



large amount of gas being photoionised by the QSO radia- 
tion. Unfortunately, spectroscopy of QSO hosts suffers invari- 
ably from the contamination of the spectrum by the emis- 
sion of the active galactic nucleus (AGN) even when ob- 
served some arcseconds away from the nucleus. Various spec- 
troscopic studies circumvente d this problem by looking at 



space-based imaging studies (e.g. 'Malkan' '1984'; Smith et aL 
1986 ; McLeod & Rieke 1994; Bahcall et a l. 1997; Dunlop et aL 

■ 120031: ISanchez & Gonzalez-SerranQll2003l) . which provide ba- 

■ sic information on the host galaxies like the morphology 
I or the luminosity. Narrow-band images centred on luminous 

emission lines have been used obtained to infer the distribu- 
tion and size of the ionised gas surrounding the QSO (e.g. 
Stockton & MacKentv 1987; Mulc haev et"ani 19961: iFalcke et al.l 



1998.:.Bennert et aL.2QQ2c.Schmitt eral.ll2003l) . However, optical 



iolJ 



and/or near-infrared spectroscopy is required for a more detailed 
investigation of the properties of the stellar population and the 

ion ised gas. 

iBoroson & Okd (Il982h were the first to detect strong B aimer 
absorption lines in the off-nuclear spectrum of the QSO 3C 48 
and highly ionised gas was alre ady found long a^o in off-nuclear 
spectra of many QSOs (e. 
119761: IBoroson & Ok3 [198" 



obscured (type 2) AGN (e.g. 'Gonzalez Delgado et al.l l200ll: 
kauff mann et al.ll2003l: ICid Fern andes et al. 2004; Ben nert et all 
20061: IStoklasova et al.l l20Q9l) 7^hich are thought to have the 
same properties as unobscured (type 1) A GN in the frame- 
work of the unificati on model of AGN (lAntonuccil 1 19931 : 
lUrry & Padovani|[l995l) . They often find a young stellar popu- 
lation in the AGN host galaxies, in particular for the most lumi- 
nous AGN, suggesting a connection between star formation and 
nuclear activity in agree ment with photometric studies of type 
1 QSO s at all redshifts (iJahnke & Wisotzki|[2QQ3l: iJahnke et al.1 
I2004allbl: ISanchez et al.ll2004bh . 



Wampl er et al. 
■Boroson et al.l 



19751: IStocktonI 
19851) indicating 



* Based on observations made with VIMOS integral field spec- 
trograph mounted on the Melipal VLT telescope at ESO-Paranal 
Observatory (programme 072.B-0550A; PI: K. Jahnke) 



iLetawe et al] (l2007l) and iJahnke et al.l (l2Q07l) . hereafter 
Let07 and Jah07 respectively, developed two different methods 
to decompose the AGN and host spectra in long-slit observa- 
tions and applied them to a sample of luminous QSOs. However, 
long-slit spectroscopy covers only a small part of the host on 
the chosen position angle and slit width. Only integral field unit 
(IFU) observations allow to study the properties of the entire host 
galaxy in an unbiased way. Dedicated AGN-host decomposition 
techniques have already been developed and were successfully 
applied to IFU data of luminous QSOs at low and high redshift 
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(ISanchez et al.ll2004at IChristensen et alJ 120061: iHusemann et alJ 
I2008h . 

In this paper, we report on optical IFU spectroscopy of the lu- 
minous radio-quiet QSO HE 1029-1401 with the VIMOS lEU at 
the Very Large Telescope (VLT). HE 1029-1401 was discovered 
by the Hamburg-ESO survey fflES: IWisotzki et al.lll99l[l2000h 
to be the brightest QSOs in the southern hemisphere with an ap- 
parent magnitude of V = 13.7 at a redshift of only z = 0.086. 

Observations with the Hubble Space Telescope re vealed a 
bright elliptical (El) host galaxy (JBahcall et al.lll997l) with an 
effective radius of r^ = 3.2'' (5.1 kpc). iJahnke et al.l (l2004ah 
found from multi-colour imaging that the host colours are sig- 
nificantly bluer than inactive galaxies at the same luminosity and 
inferred an intermediate- age stellar population of 0.7 Gyr (lumi- 
nosity weighted) by modelling the spectral Energy distribution 
(SED) with template Single Stellar Populations (SPPs). 

Spatially resolved on-nuclear longslit spectroscopy for this 
object was presented by Let07 and Jah07. Both studies consis- 
tently found highly ionised gas in the host galaxy on kpc scales. 
Jah07 found that the stellar component is non-rotating, but the 
velocity curve of the ionised gas indicated the presence of a ro- 
tating gas disc on kpc scales around the nucleus. Let07 argued 
instead that the velocity curve does not fit with a pure rotational 
velocity field based on inclination arguments. 

The main focus of this paper is to present an in-depth anal- 
ysis of the ionised gas in the host galaxy of HE 1029-1401, but 
we also perform a tentative analysis of the stellar continuum. 
Throughout the paper we assume a cosmological model with 
Ho = 70km/s, flj^ = 0-3, and Qa = 0.7. This corresponds to 
a physical scale of 1.6 kpc/'' at z = 0.086. 



2. Morphology of the host galaxy of HE 1029-1401 

A high resolution HST ima^e of HE 1029-1 401 in the F606W 
band was published by iBahcall et al.l (1 19971) . We retrieved the 
archival HST image from the Hubble Legacy Archived in order 
to re-analyse the image. An analytic point spread function (PSF) 
for the o bservations was created with the Tiny Tim software 
package (JKristl [19951) , because no sufficiently brig ht star was 
cover ed by the HST observation. We used GALFIT (JPeng et al.l 
l2002h to decompose the host and QSO c ontributions assuming a 
de Vaucouleurs light profile (Ide Vaucou leurs 1948) for the host 
plus a point source for the nucleus. The effective radius and el- 
lipticity of our best-fitting model was r^ = 3.5" (5.6 kpc) and 
e = 0.1, r espectively, which ar e consistent with the parameters 
derived bv lBahcall et"aD (Il997h . 

The nucleus- subtracted image of HE 1029-1401 is shown 
in Fig. [U A strong negative residual is still visible at the QSO 
position due to the saturated QSO core. We cleaned the im- 
age of residuals from the usual diffraction spikes by replacing 
each affected pixel with the median value of the unaffected pix- 
els within a radius of 6 pixels. Wisotzkil (1 19941) obtained red- 
shifts for 13 galaxies in the field finding 4 physical companions 
which indicate that the HE 1029-1401 is part of a loose group 
of galaxies. None of these confirmed companions were covered 
with our IFU observations. Neither CI nor C2 have redshift in- 
formation available, b ut C3 was found to be at a different redshift 
tha n HE 1029-1401 (IWisotzkil ll994. 

iBahcall et al.l (Il997 h also reported the detection of two faint 
shells located roughly 11" and 19" away from the nucleus. 
These shells are hardly visible in Fig.[T] but we highlighted their 
position by arrows for guidance and labelled them as SI and S2, 
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Fig. 1. Nucleus- subtracted HST F606W broad-band image of 
the host galaxy of HE 1029-1401. Residuals of the diffraction 
spikes have been taken out by interpolation and the image is 
smoothed by a median filter (boxsize of 4 pixels) for display pur- 
poses. Three apparent nearby companion galaxies are labelled as 
CI, C2 and C3 of which on ly C3 is confirmed to be at a diff'er- 
ent redshift (IWisotzkill 19941). Th e location of the two faint shells 
reported by .Bahcall et al.l (1 19971) are highlighted with arrows. 



respectively. Note that the shell S2 is outside the field-of-view of 
our IFU observations. 

We also performe d an isophot al analysis of the host using the 
ellipse task of IRAlQ (iTodyl 19931) after masking the 2 closest ap- 
parent companions. As the image is strongly affected by decom- 
position residuals over the central 2", we only took into account 
isophotes within a range of 2"- 10" in the semi-major axis. We 
found that the ellipticity of the host is slightly decreasing out- 
wards from 0.2-0.1 and that the a4/a Fourier coefficient is con- 
sistent with 0. Thus, the deviation from pure elliptical isophotes 
is marginal. 

3. IFU Observations and Data Reduction 

We employed the IFU mode of the VIMOS instrument 
(iLeFevre et al.ll2003l) on the VLT to perform optical 3-D spec- 
troscopy of HE 1029-1401. The observations were conducted 
on the 18th and 24th of December 2003 at a seeing of ~ 1.1", 
with the high-resolution blue and orange grisms covering the 
spectral regions around UJ3 and Ha with a spectral resolution 
of 7? ~ 2600. We chose a spatial resolution of 0.66"/fibre re- 
sulting in a field-of-view of 27"x27" that matches well with 
the angular size of the QSO host galaxy. The exposure times 
were split into 3 x 300s for the blue and 6 x 450s for the orange 
grism. A dither pattern was used in order to correct for dead fi- 
bres within the field-of-view. Additionally, a blank sky exposure 
of 300s was obtained in between of the target exposure series. 



http://hla.stsci.edu/ 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., 
National Science Foundation. 



under cooperative agreement with the 
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Arc lamp and continuum lamp exposures were acquired for each 
configuration directly after the target exposures for calibration 
purposes. Standard star exposures were taken for each night ac- 
cording to the instrumental setup. 

The VIMOS instrument is a complex IFU with 1600 operat- 
ing fibres in the high-resolution mode. These are split up into 
4 bundles of 400 fibres densely projected onto each of the 4 
spectrograph CCDs. Instrument flexure and cross-talk between 
adjacent fibres are important issues for VIMOS which need to 
be carefully taken into account in the data reduction. As this 
is not the case for the standard reduction pipeline pro vided by 
ESQ, we used our own flexible reduction package R3D ('Sanchez^ 
l2006h designed to reduce fibre-fed IFU raw data, complemented 
by custom Python scripts. The basic reduction steps performed 
with R3D include: Bias subtraction, visually checked fibre iden- 
tification, fibre tracing, spectra extraction, wavelength calibra- 
tion, fibre flatfielding and flux calibration. 

Flexure of the instrument causes substantial shifts of the fi- 
bre traces in the cross-dispersion direction between the contin- 
uum lamp exposure and the science exposure. This is a severe 
problem for an accurate extraction of the science spectra. We 
measured the fibre trace positions in the individual science ex- 
posures and continuum lamp exposure by modelling the cross- 
dispersion profile at 5 diff'erent locations on the dispersion axis 
with multiple Gaussians. The cross-dispersion profiles were gen- 
erated by co-adding 200 pixels in the dispersion direction, pref- 
erentially encompassing bright sky lines, to increase the S/N. 
The off'set between the science and continuum lamp fibre traces 
were measured to be in the range of -2.5 to 2.5 pixels. 2nd order 
Chebychev polynomials were used to extrapolate the off'sets to 
the whole dispersion axis range. 

In order to reduce the eff'ect of cross-talk we used an itera- 
tive and fast algorithm that is part of the R3D package, based 
on Gaussian profile fitting to each of the fibres in the cross- 
dispersion direction. The centroids were fixed for each fibre to 
the position on the continuum lamp taking into account the po- 
sition off'sets due to the flexure. The width of the Gaussians was 
fixed. Details of the spectra extraction algorithm can be found 
in the R3D user guidfl The width of the fibre profiles is dif- 
ferent for each of the four CCDs as it depends on the optical 
calibration of the four independent VIMOS spectrograph units. 
Thus, we estimated the fibre dispersion separately for each unit 
by simultaneously fitting each fibre with a Gaussian profile in 
the cross-dispersion direction at a single dispersion position for 
each CCD. 

The continuum lamp exposure was used to create a fibre- 
flat that corrects for the difl'erences in the fibre-to-fibre trans- 
mission. Sky subtraction and the correction for the diff'erential 
atmospheric refraction were done with custom Python scripts. 
Taking advantage of the dithered exposures we were able to 
increase the spatial sampling by a factor of 2 and to correct 
for dead fibres within the field-of-view. The flux calibrated dat- 
acubes were rescaled in flux to match with the QSO V-band 
photometry and cor rected for Galactic extinction applying the 
ICardelli et al.l (1 19891) extinc tion curve with an extinction of Ay = 
0.221 (ISchlegeletal.lll998h in the sightline to HE 1029-1401. 

4. IFU decomposition of host and QSO emission 

In order to study the emission from the stars and the gas in the 
host of HE 1029-1401 it is important to properly subtract the 
contribution of the bright QSO. Broad-band imaging studies of 
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Fig. 2. Result of the decomposition process for HE 1029-1401. 
The integrated spectrum and the decomposed QSO and host 
spectra are shown in the top panel. A zoomed version on the 
host spectrum is provided in the bottom panel. The red points 
indicate the host ma gnitudes in the V and R band measured by 
(iJahnke et al.ll2004ah . They are offset by 0.5 mag to put them on 
the same absolute scale as the IFU photometry. All spectra have 
been smoothed with a median filter of 7 pixels for display pur- 
poses. 



QSO hosts have successfully used a two-dimensional analytical 
modelling scheme to decompose the point-like nucleus and dif- 
ferent spatial extended hos t components to study the properties 
of QSO host galaxies fe.g iMcLure etal.ll 19991: iKuhlbrodtet all 



http://www.caha.es. sanchez/r3d/R3D_user_guide.pdf 



l2004l:ISanchez et al.l2004bl:lKim et al.ll2008l:lJahnke et al.l2009h . 
We have extended this method to model each monochromatic 
image of an IFU datacube in orde r to obtain the clean optical 
spectrum of the radio jet of 3C 120 (iGarcia-Lorenzo et aDl2005l: 
ISanchez et al.ll2004ai l2006l) a nd to deblend the com ponents of a 
gravitationally lensed quasar (IWisotzki et al.ll2003l) . 

Broad-band imaging studies usually estimate the point 
spread function (PSF) empirically from stacked images of stars 
within the field-of-view. The field-of-view of current IFU instru- 
ments is generally too small so that no star can generally be cap- 
tured simultaneously with the target. In light of time variable 
seeing, the QSO itself therefore needs to be used to construct 
a PSF. Fortunately, the spectral shape of the QSO is, in the ab- 
sence of atmospheric dispersion, exactly the same in each spatial 
pixel of the IFU, only scaled by a factor according to the PSF. In 
particular, the broad emission lines of type 1 QSOs are a unique 
feature of the point like nucleus, which can be used t o empiri- 
cally estimate a PSF for IFU data (iJahnke et al.ll2004cl) . 

Having obtained a PSF for the IFU data from the broad emis- 
sion lines and the morphological parameters of the host inferred 
from the HST images, we proceeded to model each monochro- 
matic IFU narrow-band image by a linear combination of the 
PSF and host model. In this way we reconstructed a pure QSO 
and host datacube containing the QSO and mean host spectrum. 
Afterwards, we subtracted the QSO datacube from the observed 
one to obtain a datacube uncontaminated by the QSO. The result 
of this decomposition process for HE 1029-1401 is presented in 
Fig.El 

We found the IFU spectrum of the QSO and host component 
to be of higher quality than the long- slit spectra presented by 
Let07. The stellar absorption lines around Mgl /15166 are clearly 
visible in our host spectrum, and the noise in the spectrum is 
greatly reduced compared to the long- slit spectrum, due to the 
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Fig. 3. Host spectrum of HE 1029-1401 in the rest- frame wave- 
length between 5100 and 5400A. The best-fitting model spec- 
trum is overplotted as the red solid line. Note that the strong 
residual of the bright [01] 5577 night sky emission line appear- 
ing at -51 35 A in the deredshifted spectrum was masked out for 
the modelling. 



much larger galaxy area captured. lEU spectroscopy improves 
the quality of the decomposition due to the large spatial coverage 
of almost the entire host galaxy and the on-source PSE estima- 
tion based on QSO features. The decomposition of the long-slit 
spectra particularly suff'ered as the only available PSE star near 
HE 1029-1401 is 2 mag fainter than the QSO (Let07). 



5. Analysis 

5.1. The host spectrum and stellar velocity dispersion 

Comparing the V an d R band apparent magnitudes reported by 
iJahnke etaP (l2004ah with the magnitudes inferred from our host 
spectrum, we found a constant off'set of 0.5 mag in both bands. 
This diff'erence can be easily explained by the brightness vari- 
ability of the QSO as we matched the overall lEU photometry to 
the archival magnitude of the QSO with a diff'erence in observa- 
tion time of several years. Note that our measured V - R colour 
in the observed frame of 0.36 mag is consistent wi th the V - R 
broa d-band colour of 0.4 ± 0.05 mag measured by iJahnke et al.l 
(l2004a) . 

While we found consistent colours with the previous broad- 
band study, we noticed a steep increase in the blue part of the 
total host spectrum. This feature is even more pronounced for 
the host spectra within 1.5'' around the nucleus. Whether this 
blue part of the spectrum originates from a young stellar popula- 
tion around the nucl eus, or from scattered QSO light due to cir- 
cumnuclear dust (e.g'Antonuc ci & Milleilll985l:IZakamska et all 
l2005l) . or whether it is an artifact of a systematic PSE mismatch 
in the decomposition process in that spectral region, we are cur- 
rently unable to say. However, we do resolve the strong absorp- 
tion lines around MglA5l66 even at the low level of signal as 
shown in Eig.O so the stellar component has a significant con- 
tribution to the continuum emission. 

We modelled the integrated host spectrum excluding the 
central 1.5'' with a linear combination of high spectral 
resolution SSP spec tra employing the STARLIGHT spec- 
tral synthesis code (ICid Eernandes et all l2005l: iMateus et al.l 
l2006h . The SSP spectra with a spectral sampling of 0.3A 



were generated with the S ed@ codqj as presented in 
iGonzalez Delgado et al.' (2005) with the following inputs: initial 
mass function from Salpeter (1955) i n the mass range .1-120 
M©; the high resolution library from Mar tins et al.l (|2005 | ) based 
on atmosphere models fro m PHOENI X (iHauschildt & BaronI 
119991: lAllard et al.l l200lh . ATLAS9 (TKuruczl Il99lh com- 
puted with SPECTRUM (iGrav & Corb allvl 1 19941) and the 
ATLA S9 library computed with TLUSTY (ILanz & Hubenvl 
l2003l) . Our input grid contains 39 SSP spectra with 
ages of 1,3,5,10,25,40,100,300,600,900,2000,5 000 and 
10 000 Myr and metallicities of 0.004, 0.02 and 0.04 Z. 

Monte Carlo simulations were used to determine the errors 
of the best-fit parameters, the flux- weighted stellar age {t^} and 
metallicity (Z) at a normalisation wavelength of Aq = 5075A, 
the extinction (Ay) and the stellar velocity dispersion (cr*). We 
generated 200 mock spectra of the host on the basis of the ob- 
served host spectrum combined with its wavelength dependent 
noise budget. These 200 artificial host spectra were analysed 
with STARLIGHT exactly as the observed host spectrum itself. 
Einally, a linear combination of only 5-7 SSP spectra was suf- 
ficient for modelling all spectra. The other SSP base spectra of 
the grid were discarded from the modelling in the first iteration 
by STARLIGHT as their total contribution to the flux density 
at the normalisation wavelength were less then 2%. The result- 
ing parameters were (h) = 5 ± 0.5 Gyr, (Z) = 0.024 ± 0.003, 
Ay = 0.5 ± 0.2 and cr, = 320 ± 90km/s (corrected for the 
instrumental resolution of ~42km/s). We repeated the same 
analysis excluding the blue part of the host spectrum (3800- 
4700A rest-frame wavelength) and arrived at consistent results 
within the errors. To our knowledge this is the first time that 
Q-^ was estimated fo r the host galaxy of such a luminous QSO. 
iJahnke et aP (l2004ah derived a younger stellar population with 
0.7-2 Gyr based on SSP modeling of the VRIJHK multi-colour 
spectral energy distribution. A direct comparison with our result 
needs to be taken with caution due to the completely diff'erent 
analysis method used to infer the characteristic age of the stel- 
lar population. However, we note that the distribution of stellar 
populations that contribute to the composite spectrum is clearly 
bimodal and consists of a young stellar population with ages 
<100 Myr (53 ± 5% of the total flux at Aq) in addition to an old 
stellar population with 10 Gyr age (47 ± 5% of the total flux at 
Aq). This indicates the presence of a young stellar component if 
the main contribution to the blue continuum of the spectrum is 
of stellar origin and not due to a featureless continuum that can 
har dly be distinguished from suc h a young stellar component 
(^e.g lStorchi-Bergmann et al.ll2000l) . 



5.2. The QSO spectrum 

The top panel of Eig. [2] shows the QSO spectrum of 
HE 1029-1401 in the full wavelength range covered by our ob- 
servations. We now focus on the analysis of the HyS-[0 iii] region 
using a multi-component fitting scheme to deblend the various 
emission lines in this spectral region. Our model consists of sev- 
eral Gaussian components and a simple linear relation to approx- 
imate the local continuum. The best-fit model to the spectrum 
and its individual components are shown in Eig.lU 

The H/S line exhibits an asymmetric line profile with en- 
hanced emission on its red wing, which can be well described by 
a redshifted Gaussian component. Erom the model we measured 



^ SedCo) is a synthesis code included in the Legacy Tool project of 
the Violent Star Formation Network; see Sed@ Reference Manual at 
ihttp : //www ■ iaa . es/~mcs/sed(a for more information 
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Fig. 4. Multi-component fit to the Hy5-[0iii] spectral region of 
the QSO spectrum. The best-fit model is indicated by the red 
solid line and the individual Gaussian components are plotted 
above the local linear continuum with the following colour cod- 
ing: Uj3 line (black), [Oiii] doublet lines (green), Feii doublet 
lines (blue). The residuals of the model are shown below with a 
refined scaling. 
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Fig. 5. Nucleus- subtracted continuum-free [O iii] A5001 narrow- 
band image. The 20 A wide band was centred on the nuclear 
[O III] A5001 wavelength. A green cross marks the position of 
the QSO. The contours indicate the position and morphology of 
the host in comparison to the ionised gas morphology. 



a line width of the broad HyS line of 6195+80 km/s (FWHM) 
and 2352+30 km/s (cr). Furthermore, the [Oiii] line profile is 
asymmetric to the blue side requiring a second Gaussian com- 
ponent. The two [Oiii] components are separated by 120+50 
km/s in rest- frame and have FWHM values 42 1+50 km/s and 
975+ 129 km/s, respectively. Additionally, a third extremely 
broad and blueshifted component is required to reach an accept- 
able fit. We are not sure whether this component is physical. 
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Fig. 6. [Oiii] channel maps created from monochromatic dat- 
acube slices corresponding to rest-frame radial velocities 
of -269.4, -180.1, -90.8, -31.2, +38.3, +87.9, +177.2, 
+266.5 km/s, respectively. The main distinct emission line re- 
gions are labelled alphabetically from A to G and the red boxes 
indicate regions that were co-added to obtain their characteristic 
spectrum (Fig. [T]). The position of the QSO is marked in each 
map by the green cross. 



representing a true extension of the asymmetric [O iii] line pro- 
file, or whether it should be rather attributed to the complex and 
asymmetric line profile of the HyS line. The continuum flux den- 
sity at 5100A (rest-frame) is (93.4+ 1 1.6) x 10"^^ ergs" ^cm'^ A" ^ 
corresponding to a luminosity of log(/lL^/[erg/s]) = 44.97 + 
0.05. 

With these measurements we estimated a black-hole mass 
from the single epoch QSO spectrum using the virial method. 
This calculation is based on the empirical relation between 
the broad-line region radius and the continuum luminosity at 
51 00 A calibrat e d bv iBentz et al.l (l2006h . and the prescription 
by ICollin et aP (l2006l) to infer the velocity of the broad line 
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Fig. 7. Co-added spectra of the regions A-G (as defined in Fig.O are shown from top to bottom split up into two wavelength ranges 
bracketing llj3 (left panels) and Ha (right panels). The black line corresponds to the observed spectrum and the red line to the fitted 
model. Two systems of Gaussians were required to model the spectrum of region G. The apparent mismatch between the model and 
^e [Nil] /16585 line in region B is due to an artifact, but the model is still constrained due to the doublet nature and the kinematical 
coupling of the lines. 



B. Husemann et al.: Mapping the ionised gas around the luminous QSO HE 1029-1401 



clouds from the UJ3 line dispersion adopting a scale factor of 
3.85. This yields a black hole mass of log(MBH/Mo) = 8.7 
for HE 1029-1401. The measurement errors for the black hole 
mass are much smaller in this case than the systematic errors 
of the method, thus we adopted a canonical error of 0.3 dex. 
Assuming the Mbr-o"* relation of iTremaine et al.l (l2002h our 
measured stellar velocity dispersion cr^ predicts a black hole 
mass of log(MBH/^o) = 8.9^^4, which is perfectly consistent 
with our virial black hole mass estimate. 

Adopting a bol ometric correction of Lboi = 9 x /1L^(5100A) 
(iKaspi et al.l l2000l) , we estimated a bolometric luminosity of 
= 46.0 + 0.2 and an Eddington ratio of 
-0.9 ± 0.2 for the QSO. 



log(LBoi/[erg/s]) 

log(LBol/^Edd) = ■ 



5.3. Ionised gas distribution and channel maps 

In Fig. [5] we present an [Oiii] /15007 (hereafter [Oiii]) narrow- 
band image (20A bandwidth) extracted from the datacube after 
subtracting the QSO contribution. Continuum emission from the 
host galaxy was previously subtracted to obtain a pure emission- 
line image. This is the first two-dimensional map of the ionised 
gas distribution around HE 1029-1401. 

The image reveals ionised gas on scales of several kpcs. 
We find a highly structured distribution of the ionised gas. The 
brightest structure is a symmetric circumnuclear bicone sugges- 
tive of an ionisation cone by the QSO. Furthermore, we detect 
ionised gas in several knots and arms and in a faint giant arc- 
like feature north-west to the galaxy centre at a projected dis- 
tance of 16kpc. Most of the distinct emission-line regions are 
well located inside the visible light distribution of the early-type 
host galaxy when compared with the HST image in Fig.[T] The 
prominent bicone-like structure is roughly aligned with the ma- 
jor axis of the host galaxy, although the orientation of the latter 
is not well constrained due to the low ellipticity of the host. 

The integrated extended [Oiii] fine flux is (328 ± 88) x 
10"^^ergs"^cm"^ corresponding to an [Oiii] luminosity of 
log(LrQjjji/[erg/s]) = 41.78 ± 0.12. This is one of the most lu- 
minous extended emission-line regions (EELRs) around a radio- 
quiet QSO c ompared to other radio-quiet, e.g. in the sample of 
IStockton & M acKenty (1987). 

In Fig. [6] we present [Oiii] channel maps corresponding to 
monochromatic images at specific radial velocities with respect 
to the rest- frame of the objects. In this way we can identify low- 
surface brightness emission-line regions that are relatively faint 
in the integrated line image or that spatially overlap but are sep- 
arated in velocity space. The main emission line regions are la- 
belled alphabetically from A to G and their characteristic spec- 
tra were obtained by co-adding several spatial pixels within the 
boundary indicated by the red boxes. These regions are a sym- 
metric bicone-like structure around the nucleus (regions C and 
G), two arm-like structures (regions D and E), two emission-line 
knots (regions B and F) and part of the giant arc-like feature 
~16kpc north-west of the nucleus (region A). We excluded the 
central -T' region around the QSO from these spectra because 
of significant residuals of the QSO subtraction. 



5.4. Spectral analysis of specific regions 

Co-added spectra of the regions A-G as defined in the previ- 
ous subsection are shown in Fig. [71 These were split up into two 
narrow wavelength ranges bracketing the UJ3 and Ha emission 
lines. No signatures of the broad emission lines of the QSO are 
visible in any of the extracted spectra after the spectral decompo- 
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Fig. 8. Diagnostic diagram for the emission-line ratios of the 
regions A-G in HE 1029-1401. The narrow emission-line ra- 
tios are indicated by black circles and the broad component 
in region G is indicated by the black opened circle. The com- 
monly u sed demarcation cur ve between AGN and Hii re- 
gions by iKewley et al.l (l200ll) is indicated as the black solid 
line . Dusty radiation p ressure-dominated photoionisation mod- 
els (iGroves et al.ll2004l) are indicated by the red tracks for four 
gas-phase metallicities with 0.25Zq, 0.5Zq, IZq and 2Zq with 
physical parameters as described in the text. 



sition. The U/3, [Oiii] ^^4959, 5007, Ua and [Nii] ^^6549, 6585 
emission lines are clearly detected in most of the selected re- 
gions. The [S 11] /1/16716, 6731 emission lines are only visible in 
the spectra of region C and G. 

We modelled each spectrum with Gaussian profiles for the 
emission lines and scaled the best-fitting stellar population spec- 
trum (cf . Sect. 15.11) to the continuum. We found that stellar ab- 
sorption lines have no significant eff'ect on the B aimer emission 
line. The flux ratio of the [Oiii] and [Nii] doublets were fixed 
to their theoretical values and all emission lines were kinemati- 
cally coupled to have the same redshift and line dispersion. We 
utilised Monte Carlo simulations to estimate realistic errors for 
all free fit parameters. 100 mock spectra were generated for each 
spectrum according to the noise of the spectrum as estimated 
from the adjacent continuum. These mock spectra were consis- 
tently analysed and the standard deviation of the resulting dis- 
tribution for each fitted parameter was taken as its Icr error. A 
3cr upper limit is given for the flux of the H/S line in region A 
as it falls below the detection limit. The measured emission line 
fluxes are given in Table [T] 

We noticed that all emission lines in the spectrum of region G 
display an asymmetric line profile. A two component Gaussian 
model for each line provided an excellent fit to the data. The ap- 
pearance of a second system of kinematically distinct emission- 
line in region G will be further analysed later in Sect. 16.41 

6. Results 

6.1. Source of ionisation and gas metaiiicity 

To infer the dominant ionisation source for the gas we em- 
ploy the commonly used emission-line ratio diagnostic dia- 
gram rOiiil ASmiimN^. rN iiU6582/Ha (lBaldwinetal.lll98ll: 
IVeilleux & Osterbrockl 19871) . The diagram distinguishes gas be- 
ing ionised by the hard UV radiation field of an AGN from gas 
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Table 1. Emission line fluxes for the UJ3, [O iii] A5001, HQf,[Nii] /16583 and the [S ii] AA6116, 6731 lines as measured from Gaussian 
fits to the spectra of the region A to G. 



Region 


H/? 


[0 III] ^5007 


Ua 


[Nii]i6583 


[S II] ^6716 


[S II] ^6731 






line fluxes in units of 10~^^ erg s~ 


icm-2 




A 


<0.5 


5.9 + 0.2 


1.9 + 0.2 


1.1+0.1 


— 


— 


B 


1.0 + 0.2 


19.5 + 0.3 


4.2 + 0.3 


4.7 + 0.2 


- 


- 


C 


7.2 + 0.2 


56.8 + 0.4 


21.7 + 0.3 


11.9 + 0.3 


4.7 + 0.3 


3.8 + 0.3 


D 


2.9 ± 0.4 


8.5 + 0.5 


4.9 + 0.5 


4.2 + 0.5 


- 


- 


E 


2.7 + 0.4 


32.0 + 0.5 


10.5 + 0.4 


6.1+0.3 


- 


- 


F 


1.3 + 0.1 


9.0 + 0.1 


3.8 + 0.2 


2.0 + 0.2 


- 


- 


G (narrow) 


5.7 + 0.5 


52.7 + 1.7 


19.1+3.4 


7.5 + 2.3 


2.8 + 1.1 


1.9 + 0.9 


G (broad) 


5.2 + 0.6 


29.9 + 1.8 


14.7 + 2.5 


10.6 + 2.2 


4.8 + 1.2 


4.5 + 1.1 



ionised by radiation from hot stars in star forming HIT regions. 
A theoretically motivated and co nservative boundary between 
both mechanisms was derived by iKewley et alJ (1200 ll) . The di- 
agnostic diagram is shown in Fig. [8] for the emission-line ratios 
inferred from the spectra of region A-G. 

We found that photoionisation by star forming regions can 
be excluded for all regions as the inferred line ratios are above 
from the Kewley et al. demarcation curve. Of course, this 
does not imply there is no current star formation in the host, 
only that AGN ionisation of the ISIM dominates strongly over 
any other ionisation mechanisms. We can also exclude shock 
ionisation based on the kinematics and the emission-line ra- 
tios. The line ratios for the broad shocked gas components in 
powerful radio galaxies are observed to h ave a much lower 
ionisation state with [O iii]/HyS~2-4 (e.g. IVillar-lMartin et al.l 
Il999|) than observed here. Furthermore, shock-io nisation mod- 
els (lDopita& Sutherland! 1 19951: lAllen et al.ir2008h indicate that 
an [O iii]/Uj3 line ratio around 10 would require shock velocities 
> 500 km/s for which we find no evidence in our observations 
(cf. Sect. 16.31) . Thus, photoionisation by the central QSO is the 
dominant ionisation mechanism for the gas in HE 1029-1401 at 
least within a radius of 16 kpc around the nucleus. 

In order to infer the metallicity of the gas we compared the 
emission-line ratios with th e dusty radiation pres sure-dominated 
photoionisation models by iGroves et aP (l2004l) . These models 
are represented by red tracks in Fig. [8] corresponding to gas- 
phase metallicities of 0.25Zq, 0.5Zq, IZq and 2Zq assuming 
an electron density of ^ = lOOcm"^ for the emitting clouds and 
a power-law index a = -1.4 for the ionising QSO continuum. 
The line ratios do not strongly depend on the electron density, 
so we chose the model grid with the lowest density close to the 
estimated value as discussed in the next section. Note that the 
true metallicities of the photoionisation models are a factor of 2 
higher than the gas-phase metallicities, since half of the overall 
metal content is depleted into dust in the models. The compari- 
son with the photoionisation models revealed that the emission- 
line ratios follow quite closely the model with a solar metallicity. 
Only the gas in region B appears to have a systematically diff'er- 
ent ionisation state and/or metallicity. 



6.2. Total ionised gas mass 



To estimate the total ionised gas mass in HE 1029-1401 we fol- 
lowed the prescription and assumptions of lOsterbrock & Ferlandl 
(l2006l) . By measuring the luminosity of the H/S recombination 



line one can estimate the ionised gas mass Mic 
lowing formula: 



Mion = 



1.4mn 



Hea^hVufi 



^Uj3 ■ 



using the fol- 



(1) 



where rup is the proton mass, h is Planck's constant. He is the 
electron density, a^ is the eff'ective recombination coefficient 
for UJ3 and L^s is the HJS line luminosity. Assuming Case B re- 
combination, the low-density limit and a typical gas temperature 
of r ^ 10, 000 K for ionised nebulae, we end up with 



Mi, 



10^ 



100 cm" 



--Uj3 



1041 erg s-1 



(2) 



In our case we were able to measure the electron density 
sensitive [Sii] /16716//16731 emission-line ratio only in high- 
surface brightness regions of the bicone structure (region C 
and G). The line ratios in region G are rather uncertain due 
to the line blending of a narrow and broad component, so 
that we estimated the electron density from region C only. We 
measured a [Sii] /16716//16731 line ratio of 1.2 ± 0.1 for re- 
gion C. This ratio implies an electron densities of 200^|^Qcm"^ 
(lOsterbrock & Fer land 2006, p. 123) assuming an electron tem- 
perature of 10 000 K. 

Since the narrow Hfi line is the weakest in all of the spectra, 
we used the [O iii] emission line as a surrogate for the UJ3 emis- 
sion assuming a line ratio of [Oiii]/Ufi ~ 10 (cf. Fig. (8]). With 
the integrated [Oiii] line luminosity measured from the [Oiii] 
narrow-band image (cf. Sect. l5.3b and our adopted electron den- 
sity, we estimated an ionised gas mass of Mion = S'^^^IO^ M©. 
Note that we did not take reddening due to dust in the host galaxy 
into account, which would increase Lj^ and hence the ionised 
gas mass. Our inferred ionised gas mass should be taken as a 
rough estimate due to the various assumptions needed and can 
only be used at an order-of-magnitude level. 



6.3. Global gas kinematics 

In order to obtain full 2-D velocity and velocity dispersion maps 
of the ionised gas we fitted the [O iii] doublet lines as they are 
the brightest emission lines in all of the individual 6400 spectra 
of the nucleus- subtracted datacube. The measured emission-line 
centroids were converted into radial velocities with respect to 
the systemic redshift for each spatial pixel to construct a 2-D 
velocity map (Fig. [9] left panel). A different visualisation of the 
velocity field are 'long-slit' velocity curves (Fig. [9] right panel) 
extracted from our velocity map for the position angles (PAs) 
25° and 60°. We chose these two PAs such that one shows the 
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Fig. 9. Left panel: [Oiii] velocity map of 
HE 1029-1401 with respect to its estimated 
systemic velocity (25670 km/s). The two black 
lines correspond to position angles of 25° 
and 60°, respectively. Right panel: Velocity 
curves extracted for these two orientations 
from the maps. Negative offsets correspond to 
the South-East side and positive offsets to the 
North- West side of the nucleus. 
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Fig. 10. Same as Fig. [51 but for the [O iii] ve- 
locity dispersion, corrected for the spectral res- 
olution of the instrument. 



velocity curve of the bi-cone structure and the other to cross 
the arm-like feature. Similarly, we present a 2-D [O iii] veloc- 
ity dispersion map (Fig. [TOl left panel), corrected for the spectral 
resolution of the instrument, and the corresponding 'long- slit' 
curves (Fig. [TOl right panel). The gaps in the velocity and dis- 
persion curves around the circumnuclear region are due to the 
strong residuals of the decomposition process which prevent a 
reliable measurement of the extended emission lines. 

The 2-D kinematic information from our IFU data draws a 
much clearer picture than previous long- slit studies of this ob- 
ject. While the [O iii] velocity curve at a PA of 25° appears like 
a symmetric rotation curve up to 5'' from the galaxy centre, it is 
highly disturbed at larger distances. At a PA of 60° the velocity 
curve crosses the arm-like structure detected in the [O iii] flux 
distribution which clearly breaks the symmetry seen at a PA of 
25°. We measured a maximum amplitude of ±170 km/s for the 
radial velocities at a distance of 4.8 kpc (3'') from the nucleus, 
which is consistent with the measurements by Let07 and JahOT 
based on long-slit spectroscopy only. However, the long-slit ve- 
locity curve presented by Let07 covers only ±4'' which is less 
than half of the region we studied with our IFU data. Therefore, 
they miss the turnover of the velocity curve at 4''-6''from the 
nucleus. Another major kinematic component is the prominent 
distortion in the rotational velocity field on the West side of the 
nucleus, which was missed by the long-slit data of LetOT and 
Jah07. 

The measured velocity dispersion {cr) of the gas after cor- 
recting for the instrumental resolution is ~ 30 km/s on the red- 
and ~40 km/s on the blue- shifted side of the velocity curve. On 
both sides the vjcr ratio is significantly greater than 1, which in- 
dicates rotational support of the gas in a disc. This is consistent 
with the conclusion of JahOT. The velocity dispersion increases 
significantly up to ~ 80 km/s between the region of the gas disc 
and the major kinematic distortion on the West side. This is most 
probably related to the superposition of the two kinematic com- 
ponents as we find asymmetric and double-peaked emission- 
lines in this particular region. In the circumnuclear region the ve- 
locity dispersion rises steeply towards the centre. We will study 
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Fig. 11. Detailed modelling of the emission lines in the spectrum 
of region G. The emission lines around HyS are shown in the top 
panel and the ones around Ha in the bottom panel. The spectrum 
and its model together with the individual Gaussian components 
for the emission lines are plotted with colours given the in leg- 
end. 



the kinematics of this region in greater detail in the next section. 
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Fig. 12. Distribution and kinematics of two kinematically distinct emission line components. The flux distribution (left), radial 
velocity (middle) and velocity dispersion (right) are shown for the narrow emission line in the top panels and in the bottom panels 
for the broad component. 



6.4. Gaseous kinematics in the inner3kpc 



As previously mentioned in Sect. 15.41 a two component model 
for the emission lines was required to model the spectrum of 
region G. We present the model for that spectrum and the indi- 
vidual emission-line components in Fig.[TT] All the strong emis- 
sion lines are clearly asymmetric and a second broader Gaussian 
component is needed to match this asymmetry. Furthermore, we 
found that all emission lines share the same kinematics and were 
therefore kinematically coupled in the model. We inferred a ve- 
locity dispersion of 45 km/s for the narrow and 110 km/s for the 
broad component. 

From the spatially resolved kinematics assuming only a sin- 
gle Gaussian for the emission line, we note that the velocity dis- 
persion is rising close to the centre. This suggests that the broad 
component is also spatially resolved. Thus, we modelled each 
spectrum of the datacube with a two-component Gaussian model 
for the [O iii] doublet emission line. We employed the statistical 
F-test on the x^ values for the two models to decide weather the 
two component model is a statistically significant better repre- 
sentation of the spectrum than a single component model. The 
resulting maps of the flux distribution, radial velocity and veloc- 
ity dispersion for the narrow and broad components are shown 
in Fig. [H 

These maps reveal that the narrow component originates 
from rotationally supported gas without a major increase in the 
velocity dispersion as close as 1 kpc to the QSO. The broad com- 
ponent, on the other hand, is almost non-rotating with a velocity 
gradient of only ±40 km/s compared to ±120 km/s of the rotat- 
ing gas in the same region. The spatially resolved velocity dis- 
persion of the gas in the broad component were measured to 



be in the range of (Xg = 100 - 140 km/s in agreement with the 
velocity dispersion inferred from the integrated spectrum of re- 
gion G. Because cr/v is greater than 1 for the broad component, 
this gas seems to be dispersion dominated and mainly supported 
by random motion. We also detected the broad component on 
the blueshifted side of the rotating gas coincident with region C, 
which is hardly noticeable in the integrated spectrum of that re- 
gion as the flux in the broad component is much lower than flux 
in region G. This dispersion dominated gas component in the 
circumnuclear region was not detectable in the previous studies 
by Let07 and Jah07 due to the lower spectral resolution. 



6.5. Surface brightness profiies and apparent neighbouring 
galaxies 

We measured the radial surface brightness profiles of the stel- 
lar continuum emission as seen by HST and of the [O iii] line 
emission from the VIMOS data using circular annuli. To ac- 
count for the biconical structure of the [O iii] light distribution 
we also computed surface brightness profiles within a 70° wide 
cone along the major and minor axis. Due to the low spatial res- 
olution of the VIMOS date we took the subpixel coverage of the 
annuli and cones into account. Since the ratio of [O iii] to HJS is 
constant of -10 for the entire galaxy, it is reasonable to assume 
that the [O iii] light distribution is a tracer of the ionised gas dis- 
tribution. The profiles are compared in Fig. [13] after the HST 
surface brightness was scaled to match the total [O iii] surface 
brightness in the radial distance range from 7.5'' to 9'\ 

Surprisingly, the surface brightness profile of the total 
ionised gas apparently follows that of the stars on a global scale 
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Fig. 13. Radial surface brightness profiles of the stellar contin- 
uum (solid line) compared with the total ionised gas distribu- 
tion (dashed line), the ionised gas within a 70° wide cone along 
the major axis of [O iii] (dotted line) and the ionised gas for a 
similar cone along the minor axis (dashed-dotted line). The sur- 
face brightness profile of the stellar continuum was scaled to 
match that of the total ionised gas between 7.5'' (12 kpc) and 9'' 
(14.5 kpc) radial distance. The radial distances of the apparent 
companions CI and C2 are marked for comparison. 



despite the morphological diff'erence. However, the profiles are 
significantly diff'erent for the light distribution in the cones along 
the major and minor axis. The drop in the [O iii] surface bright- 
ness along the minor axis is steep with a power-law index of - 1 .8. 
This is roughly consistent with the geometrical decrease in the 
ionising photon density that scales as r"^. On the other hand, the 
profile along the major axis is much shallower also with respect 
to the stars indicating that the radial ionised gas distribution is 
determined by the gas density rather than the amount of available 
ionising photons in this direction. The distinct [Oiii] emission 
line knots and arcs as identified in the [O iii] narrow-band image 
and channel maps (cf. Fig.[5]and Fig. [6]) produce excesses in the 
total [O III] surface brightness profile and along the major axis 
at around 5 kpc, 10 kpc and 16 kpc radial distance as expected. 
We do not detect any corresponding features in the continuum 
surface brightness profile. The stellar surface brightness profile 
itself is rather smooth except for the bright apparent companion 
CI located 4arcsec away from the host centre, but we do not 
find any significant excess in the [O iii] surface brightness at that 
location neither in the major nor the minor axis cone. 

We consider two explanations for the lack of ionised gas at 
the location of the bright apparent neighbouring galaxy CI north 
to the nucleus: 1.) The galaxy is not a physical companion so its 
gas cannot be ionised by the QSO. 2.) The galaxy is a physical 
companion with an intrinsically low amount of gas to be ionised. 

The relation of the second faint galaxy to the host galaxy of 
HE 1029-1401 remains unclear. We detected weak [Oiii] emis- 
sion at its position which corresponds to one end of the large 
arc-like structure seen at a distance of 10''. It is possible that 
gas is stripped off' from the faint companion while it is orbiting 
around the host galaxy forming a tail of gas that is partially il- 
luminated by the QSO radiation. Since redshift information is 
not available for this galaxy, a simple projection eff'ect is also a 
reasonable explanation. 



7. Discussion 

7.1. The unexpected anisotropic AGN radiation field 

The biconical distribution of the ionised gas in the circumnu- 
clear region of HE 1029-1401 appears to be very similar to 
the ionis ation co nes detected i n several nearby Seyfert galax- 
ies (e.g. IPoggd [l9 88a.b: Stor chi-Bergmann & Bonattol 1 19911: 
IStorchi-Bergmann et al.lll992l) . These collimated AGN radiation 
cones are supporting the unified model for AGN in which the 
nucleus is thought to be surrounded by a dusty torus providing a 
partial obscuration of the AGN. This model explains the distinc- 
tion between type 1 and type 2 AGN solely by the orientation of 
the obscuring torus with respect to the observers line-of- sight. In 
this picture, the opening angle of the ionisation cones is simply 
determined by the covering factor of the torus. 

Why is such an anisotropic radiation field on kpc scales 
unexpected for a luminous type 1 QSO like HE 1029-1401? 
Firstly, ionisation cones are best seen if they are oriented per- 
pendicular to our line-of- sight, which means that the torus is 
oriented almost edge-on and the AGN is obscured. This is 
the reason why most of the known ionisation cones were de- 
tected in Seyfert 2 galaxies. Secondly, the 'receding torus' model 
(iLawrencdl 19911) predicts the opening angle of the AGN ionisa- 
tion cones to increase with the AGN luminosity, which is ob- 
servationally supported by the decreasing fraction of obscured 
AGN wit h increasing luminosity in X-ray selected AGN sam- 
ples (e.g. IUedaetal.1 [20031: iHasinged [20041 l2008h . The X-ray 
2-10 keV luminosity o f HE 1029-1401 is 2.3 x lO^^ergs'^ 
(^Re eves & Turnej2000l converted to a concord ance cos mology) 
corresponding to a type 2 fraction of ^30% (iHasinger 2008) . 
From this fraction one would expect a large an opening angle of 
-115°, much more than what our observations suggest. 

iMulchaey et all (Il996l) simulated the apparent emission-line 
images of the ionisation cones assuming that the ambient gas dis- 
tribution is a spheroid or disc. They found that for a spheroidal 
gas distribution the emission-line image may have a V-shape 
morphology but the apparent opening angle is always equal or 
greater than the opening angle of the ionisation cones. Due to 
the rotational motion of the gas, our observations indicate that 
the gas is more likely to be distributed in a disc rather than 
in a sphere. In that case, the apparent opening angle from the 
emission-line images is always smaller than the opening angle 
of the ionisation cones depending on the angle between the disc 
and cone axis. If the ionisation cones is oriented such that it in- 
tercepts the disc almost at its edge, the apparent opening angle 
would be much smaller and consistent with our observations. 
Thus, the morphology of the emission-line gas around HE 1029- 
1401 can be consistent with the expectation of the unifed model 
of AGN if the ambient gas is distributed in a disc combined with 
a certain configuration of the AGN ionisation cones with respect 
to the disc. 



7.2. The origin of the gas 

We inferred from our IFU data that the elliptical host of 
HE 1029-1401 contains ionised gas with a mass of the order 
of 5 X 10^. The ionised gas mass is rather small compared to 
the stellar mass of the system that has been measured to be 
^stellar ~ 1 X lO^^M© bascd on broad-band SED modelling (M. 
Schramm, private communication). Unfortunately, the neutral 
gas fraction of HE 1029-1401 is unknown. Thus, it is unclear 
if HE 1029-1401 contains more gas than comparable inactive 
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early-type galaxies, or whether the luminous QSO is just able to 
ionise a significant fraction of the neutral gas content. 

But where did the gas come from? Our emission-line diag- 
nostic analysis showed that the metallicity of the gas is close to 
solar. We now compare this with the mass- metallicity relatio n of 
star- forming galaxies as presented by Tremonti et al.l (|2004|) . As 
the mass-metallicity r elation depends criti cally on the assumed 
metall icity estimator (JKewley et al.l l2006l) . the iTremonti et al.l 
(|2004|) relation best matches with our metallicity esti mation as 
it is a lso based on the CLOUDY photoionisation code (JFerlandl 
Il996l) . What we find is that the expected value for the oxy- 
gen abundance for a galaxy with a stellar mass of 10^ ^Mq is 
12 -\- log(0/H) = 9.1 1. This is more than 0.3 dex larger than the 
solar oxygen abundance measured for HE 1029-1409. Although 
the scatter in the mass-metallicity relation is non-negligible, the 
metallicity of the gas in HE 1029-140 9 woul d correspond to a 
>3(T deviation from the lTremonti et al.l (|2004|) mass-metallicity 
relation at a stellar mass of IO^^Mq. 

iPeeples et aJ] (l2009h and lAlonso et al] (l2010h recently 
showed that the high-mass low-metallicity outliners of the mass- 
metallicity relation are morphologically disturbed galaxies with 
bluer colours (or younger stellar populations) suggesting an in- 
flow of low-metallicity gas driven by galaxy interactions and 
accompanied by enhanced star formation. This is consistent 
with the observations of HE 1029-1401, but the signatures for 
ga laxy interactions are weak with only the faint shells detected 
by lBahcallet al.l (1 19971) . The high stellar mass and smooth sur- 
face brightness profile of the host suggest that the last major 
merger that formed the elliptical galaxy happened at a much 
earlier epoch. Assuming that the gas originates entirely from a 
single progenitor galaxy, the solar gas phase metallicity would 
correspond to a stellar mass of roughly 5 x IO^Mq almost 10 
times smaller than the host mass of HE 1029-1401 itself. This is 
consistent with the picture that ionised gas originates from one 
or a few minor companions that merged with HE 1029-1401. 
Therefore, we favour to interpret the low metallicity of the 
ionised gas in HE 1029-1401 as a clear indication for an ex- 
ternal origin of the gas via the infall of one or a few minor com- 
panions. 

Similar results have been obtained for the EEL Rs around 
the radio-lou d QSOs 4C 37.43 (iFu & Stocktonll2007h and 3C 79 
(IFu & Stockton.2008.) . The authors found in both of these cases 
a solar to sub- solar metallicity of the ionised gas on kpc scales 
around their massive elliptical host galaxies containing black 
holes with a mass of the order of lO^M©. They also came to 
the conclusion that the metal-poor gas needs to have an external 
origin. 



7.3. The kinematics of the ionised gas - Indications for 
mergers or AGN outflow? 

The [O III] emission-line profile in the nuclear spectrum is asym- 
metric, which we resolved into two distinct Gaussian com- 
ponents with a FWHM of 421±50km/s and 975±129km/s, 
respectively, where the broader component is blueshifted by 
120 ± 50km/s. Such asymmetric [Oiii] lines have been pre- 
vious l y reported in many oth er QSOs (e.g. iHeckman et all 
1 19811: IVeron-Cetty et al.l l200ll) and were interpreted as an 
indication for high ionisa tion outflows from the nucleus 
(IZamanov & Marzianill2002h . We can not spatially resolve the 
extended ionised gas closer than ~1 kpc from the nucleus, but 
outside this unresolved region we did not find any signatures for 
a similarly broad and blueshifted ionised gas component indi- 



cating that this possibly outflowing component is restricted to 
scales below ~1 kpc. 

Instead, we find that the majority of the extended ionised gas 
appears to be bound in a rotationally supported gas disc. The 
lack of rotational motion in the corresponding stellar component 
as inferred by Jah07 further supports an external origin of the gas 
from a kinematically point of view. The non-rotating dispersion- 
dominated ionised gas within the central 3 kpc could either rep- 
resent an intrinsic reservoir of gas originating from the stars in 
the host or gas that is directly heated by the AGN. 

On larger scales the kinematics of the ionised gas around 
HE 102 9-1401 strongly de viates from the rotational motion of 
the disc. lGarcia-Lorenzo et al. (2005) found a similarly distorted 
kinematic pattern in the gaseous velocity field and inner shells 
around the nucleus of 3C 120, which they interpreted as evi- 
dence for a merging event. Above we invoked minor mergers 
to explain the external origin of the gas, which would also pro- 
vide a consistent explanation for the distortion in the velocity 
field. On the other hand, morphological distortion in the stellar 
component should be detected at the location of the kinemati- 
cal distortions if mergers play the dominant role. It could be that 
the current HST image of HE 1029-1401 is simply too shallow 
to detect faint morphological distortions. iBennert et al.l (l2008l) 
found shells and tidal tails in four out of five elliptical QSO 
host galaxies observed with HST in deep exposures suggesting 
that merger are quite common for such galaxies. But without a 
matched control sample of inactive elliptical galaxies this result 
cannot be taken as a direct evidence for a causal link between 
mergers and QSO activity. Indeed, in a large sample o f mainly 
inactive elliptical galaxies presented by lTal et aP (l2009l) , similar 
distortions were found in 75% of them with no correlation be- 
tween AGN activity and large scale distortions. However, their 
AGN comparison sample for comparison consists mainly of low- 
luminosity radio galaxies which might have a diff'erent fuelling 
mechanisms compared to luminous QSOs. 

As an alternative possibility we speculate that these struc- 
tures could be ionised filaments on the surface of an expand- 
ing shell/bubble driven out by the AGN. Two examples of 
such outflowing bubbles were investigated recently with de- 
tailed IF U observations of the Broad Absorption Line QSOs 
Mrk 231 (lLiparietal.1 120091) and HE 0450-2948 (iLipari et al.l 
l2009l) . They argued that these outflows are part of galactic out- 
flows driven by extreme starbursts as the result of major mergers 
which possibly lead also to the onset of QSO activity due to an 
interaction between the nuclear star formation and the central 
black hole. Such a particular phase of QSO activity has been 
claimed to be part o f an evolutionary sequence of AGN (e.g . 
Sanders etal."1988a''b'; 'Lipari"1994'; 'Canalizo & Stockton|l200d; 
Papadopoulos et al. 2008; Lipari & Terlevich 2006, and refer- 
ences therein) as a short transition phase before the QSO be- 
comes dominant. In this scenario HE 1029-1401 would have to 
be place at the end of a potential AGN evolutionary path, the 
QSO dominant phase. Yet, due to the unknown geometry and 
proper motion of these intriguing structures, it is difficult to pin 
down the nature of the disturbed velocity field as due to an ex- 
panding bubble. An interpretation as signatures of (minor) merg- 
ing appears at least as likely. 



8. Summary and Conclusions 

We presented optical integral field spectroscopy observations of 
the luminous radio-quiet QSO HE 1029-1401 with VIMOS in- 
strument at the VLT. A dedicated decomposition technique for 
IFU data was used to decompose the spectra of the host and 
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the QSO component. This allowed us to perform a detailed 
spatially-resolved analysis of the ionised gas properties based 
on emission-line diagnostics and a tentative analysis of the in- 
tegrated stellar continuum at a S/N level of ~10. We draw the 
following conclusions from our present study of this QSO: 

- We found tentative evidence from the stellar continuum for 
a blue featureless continuum or a young stellar population 
(<100Myr) embedded in the old stellar population (10 Gyr) 
of the massive elliptical host galaxy. The stellar velocity 
dispersion is estimated to be cr* = 320km/s, although the 
uncertainty of ±90 km/s is considerable. 

- The black hole mass is estimated from QSO spectrum to 
be log(MBH) = 8.7 + 0.3 accreting at an Eddington ratio of 
log(LBoi/^Edd) = -0-9 ± 0.2. Our Mbh and cr* measurements 
are co nsistent with the Mbh - cr* relation of iTremaine et al.1 
(l2002h within the uncertainties. 

- We detected highly ionised gas in the host galaxy up to a 
projected distance of at least ~ 16kpc that is solely ionised 
by the QSO radiation. The integrated [Oiii] luminosity of 
this EELR is with log(L([Oiii])/[erg/s]) = 41.78+0.1 one of 
the most luminous EELR detect around a radio-quiet QSO. 
The total ionised gas mass was estimated to 3^2 >^ IO^Mq, 
which should be correct at an order-of-magnitude level only. 

- The morphology of the ionised gas appears to be biconical 
with additional faint knots and arcs at larger distances. The 
rather narrow bicones may imply more collimated ionisation 
cones of the AGN than predicted by the unification model 
for such a luminous QSO. Assuming a disc-like distribution 
of gas, we prefer to interpret the biconical emission-line 
image as the projection of the partially illuminated gas disc. 

- Comparison of photoionisation models with the observed 
emission-line ratios of the gas indicated that the metallicity 
of the gas is significantly lower than expected from the 
mass-metallicity relation of galaxies. We interpret this as an 
external origin of the gas and argue that this is most likely 
due to minor merger(s). 

- We found that the ionised gas is rotationally- supported 
with an additional dispersion-dominated component in the 
circumnuclear region up to ~3kpc. The ordered velocity 
field is clearly distorted on larger scales that either support 
the scenario of minor mergers of the gas or might be 
speculatively interpreted as filaments on the surface of an 
AGN-driven outflowing shell/bubble. However, we did not 
find clear and conclusive evidence for an AGN outflow. 



Our observations draw a scenario of HE 1029-1401 in 
which the massive elliptical host has recently accreted fresh gas 
from its environment via minor mergers. This process leaves the 
host morphology almost unchanged, but the black hole accretion 
can be 're-juvenated' and a disc of gas is formed. If the accre- 
tion event is accompanied by a short burst of star formation, the 
previous red elliptical galaxy would become bluer moving the 
galaxy back to the 'green valley' in the colour-magnitude dia- 
gram. The importance of minor mergers in the late evolution of 
the early-type galaxies population has re cently been high lighted 
(iNaab et al.fl2"009l: [Johansson ^"10120091: iBezanson et al] [2009). 
It might also be an important process to explain a luminous AGN 



phase in early-type galaxies as indicated by our single-object 
study. 
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